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Abstract 

Searching and sorting used as a subroutine in many important algo- 
rithms. Quantum algorithm can find a target item in a database faster 
than any classical algorithm. One can trade accuracy for speed and find a 
part of the database (a block) containing the target item even faster, this 
is partial search. An example is the following: exact address of the target 
item is given by a sequence of many bits, but we need to know only some 
of them. More generally partial search considers the following problem: a 
database is separated into several blocks. We want to find a block with 
the target item, not the target item itself. In this paper we reformulate 
quantum partial search algorithm in terms of group theory. 

1 Introduction 

Database search is used as a subroutine in may important algorithms |16[I14|[T^ . 
Grover discovered a quantum algorithm that searches faster than any classical 
algorithm ^ . It consists of repetition of the Grover iteration Gi . We shall call 
it global iteration, see lO. The number of repetitions [queries to the oracle] is: 

Jfull = f^ as N^oo (1) 

for a database with N entries. There is no faster quantum algorithm 01|S1[7]- 
Nevertheless if we need less information then the search can be performed faster. 
For example if the exact address of the target item in a database is given by 
a sequence of many bits 6162 •■• but we want to know only three of them, 
we can do it faster then iQJ. This is an example of partial search. A partial 
search considers the following problem: a database is separated into K blocks 
of a size h = N/K . We want to find a block with the target item, not the 
target item itself. The block with the target item is called the target block, 
all other blocks are non-target blocks. Grover and Radhakrishnan suggested a 
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quantum algorithm for a partial search in [S]- Partial search naturally arise in 
list matching Classical partial search takes ~ {N — b) queries, but quan- 
tum algorithm takes only ~ {^/N — cocSVb) queries. Here coefF is a positive 
number, which has a limit, when number of items in each block is very large 
6 — > oo. Grover-Radhakrishnan algorithm uses several global iterations G^^ and 
then several local iteration , see ((HJ ■ Grover-Radhakrishnan algorithm was 
simplified and clarified in and optimized in [the number of queries to 
the oracle was minimized, positive coefF was increased]. Other partial search 
algorithms were studied in The algorithm for blocks of finite size was for- 
mulated in [S]. Partial search algorithm for a database with multiple target 
items was formulated in |13| . 



2 The GRK Algorithm for Partial Search 

To introduce the partial search algorithm it is useful to first remind the full 
Grover search. We shall consider a database with one target item. The Grover 
algorithm finds the target state \t) among an unordered set of A'' states, which is 
called the database. In the classical case, the items in the database are labeled 
by a sequence of bits. In the quantum case these sequences label orthonormal 
basis in a linear space \x). The search is performed by repeating global iteration 
which is defined in terms of two operators. The first changes the sign of the 
target state \t) only: 

it^i-2\t)(t\, m = i, (2) 

where / is the identity operator and \t){t\ projects on the target item. The 
second operator, 

4 =/-2|si)(si|, (3) 
changes the sign of the uniform superposition of all basis states |si) 



N-l 
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{s^\s^) = ^. (4) 



x=0 

The global iteration is defined as a unitary operator 

Gi^^isjf (5) 

To describe partial search we need to introduce a local search. Local search 
is a search inside of each block done simultaneously in all blocks. In one block 
local search acts as 

G2 = -isjt, (6) 

where It is the same operator appearing in the global search (0), but Ig^ acts 
in each block as 



1 

one block 
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Local search in the whole database is the direct sum of /s^ with respect to all 
blocks. 

Let us mention eigenvectors and eigenvalues of global iteration, see 



\ 



(8) 
(9) 



where 9i is 



sin^ 6*1 



In the case of local search G2 they are given by 



-^|i)±-^|ntt) (10) 



Amplitudes of all non-target items in the target blocks are the same. So we can 
only follow the amplitude of |ntt) which is the normalized sum of all non-target 
items in the target block: 



ntt = 



(nttlntt) = 1, (nttlt) = 0. 



(11) 



target block 



Here the angle 62 is given by 



.2. K 1 
sm 92 = — = - • 
N b 



(12) 



2.1 Three Dimensional Space 

Amplitudes of all items in non-target blocks are the same. So we can only follow 
the amplitude of \u). The vector \u) is given by 
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E 



(uu) = 1. 



(13) 



all items in all 
non-target blocks 



Together with the target item \t) and \ntt) [sec (jllf) ] the unite vector \u) form 
a orthonormal basis in three dimensional linear space: {t\u) = {ntt\u) — 0. 
For example the eigenvectors of global iterations Gi in Eq. (8) can be written 

as 



(14) 
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Below we shall describe the partial search algorithm as a three dimensional 
matrix if^ . 

Meanwhile let us remind the GRK algorithm for partial search. The partial 
search of JTT1 creates a vector 

\d)=GiGi'G{'\s,). (15) 

Note that this algorithm uses a sequence of global- local operators. The final 
operation Gi is necessary since G2 acts trivially on \u), i.e. {u\G2 = {u\. The 
final state \d) should have zero amplitudes of each item in non-target blocks, in 
other words it should satisfy 

{u\d) = 0, 

which means that a measurement will reveal the position of the target block. 
We consider large blocks b = N/K 00. The number of blocks K is a. finite 
number {K — 2,3, ... ), it is an important parameter in the algorithm. It is 
useful to introduce an angle 7 defined by 



In the limit 6 — > 00 it was shown in ^ that the total number of iterations 
grows as 



where the coefficients rjK and uk have a well defined limit. The total number 
of queries to the oracle is ji -I- ^2 + 1 — > tt/A^/N — {rjK — aK)\/N/K. 

One of the authors found the optimal values of (riK, ax) (see such that 
total number of queries to the oracle is minimal: 



tan —= ~ , cos[2aK) ^ 



kJ K-2 ' ^ ' 2{K-1) 



3 0(3) Group Formulation of GRK 

The GRK is the fastest among partial search algorithms, which use the sequence 
GiG2^&i. Partial search algorithms using some other sequences were consid- 
ered in but no acceleration was found. In order to prove that the GRK 
algorithm is the fastest among all partial search algorithms, which use other 
sequences of Gi and G2 one has to prove the following: 

Conjecture: Start from an arbitrary vector |0) in the three dimensional space 
with the basis (|t), \ntt), \u)). If the sequence of local-global-local-global opera- 
tions can find the target block 

(zi|GfG?GfG^2M0> =0, 
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then there exists a global-local-global sequence such that 

(u|GfG^^Gi^|0) = O, 
it finds the target black faster, i.e. + j2 + ji < ji + j2 + js + ji- ^ 

If the conjecture is true then GRK is the fastest among partial search algo- 
rithms using arbitrary sequence of local and global searches. 

Proof: Assume a general algorithm constructed with global Gi and local 
G2 iterations with n letters given by 

(M|Gf G^2"-^Gf -^G^2"-'Gf -"G^2"-' • • • Gf G^2'Gf |si) = 0, (16) 

where ji > for i < n — 1 and j„ > 0. In particular we can have ji — 0, so the 
sequence can also start with G2. The final operation has to be Gj" since G2 
acts trivially on \u), i.e. (it|G2 = {u\. Equation l|16|) can be written as 

where = Gj""" • • • G2^G{^|si) Using the conjecture, this can be reduced to 

(u|Gf-Gh^Gf-^|0) = O, 

where j„_i -I- jn-2 + jn-3 < jn + jn-i + jn-2 + Jn-4- We uow have another 
sequence 

(u|Gf-G^"-Gf-Gf-G^"-|0') = 0, 
where |(/)) — G{""''G2""^ We can put two powers of Gi together 

{u\Gl-^Gt^Gi-^+'-^&r'^W) ^0, 

and use the conjecture again, starting with After several iterations, H16|l 
will be reduced to 

(w|G^G^G^^G^|si) =0. (17) 

It was proved in ^21 that the GRK algorithm is the fastest among all pos- 
sible algorithms in the form H17|) . This means that in order to proof that GRK 
algorithm is the fastest among all partial search algorithms consisting of arbi- 
trary sequence of of Gi and G2 it is enough to prove the conjecture. ^ 

• Let us reduce partial search to 0(3) group. 

In subsection 2.1 we explained that the partial search algorithm acts nat- 
urally in three dimensional space with the orthonormal basis; target item \t), 
non-target items in the target block \ntt) and all items in non-target blocks \u). 
Search operations are rotations in three dimensional space spanned by these 
three vectors. All the vectors involved in present quantum search problem can 
be written in this basis as 

\V)^\ 6 1 , (18) 
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In the above equation (a, b, c) are the real coefficients in the base defined by 
\ntt), \u)), meaning 

\V) = a\t) + h\ntt) + c\u) (19) 
For example, the initial state Q can be written as: 

sin 7 sin 62 

|si) = I sin7COs6'2 I (20) 



cos 7 



and the local uniform state {T)) is 



sm 62 

\S2)=\ cos ^2 I (21) 

From this basic relations and the definitions of G2 we can calculate its j2 
power of local search: 

/ cos(2j2^^2) sin(2j202) \ 
G{'=\ -sin(2j2e2) cos(2j2e2) (22) 
V l) 

The ordering of eigenvectors is \ntt) and \u). The matrix has three eigen- 
vectors: 

\v^) = eM±'^^02j2)\vt), G? \vl) = K) (23) 
Where the eigenvectors can be written as 



L j , \vl) = f I • (24) 



V2 

In the same way, ji repetitions of the global iterations lO is 

= (25) 



c(2ji0i), s(2ji0i)s(7), s(2ii0i)c(7) 

-s(2ji0i)s(7), (-l)^^c2(7) + c(2ji0i)s2(^), s(7)c(7) ((-1)^^+' + c(2ji0i)) 
-s(2ji0i)c(7), s(7)c(7) + c(2ji0i)) , (-l)^^s2(7) + c{2ne^)c\^) 

This is a simplified asymptotic expression valid in the limit of large blocks 
6 — > cxD. The matrix has three eigenvectors 

Gi>fHexp(±2*0iji)K±), G\'\vl) = {-ir\vl) (26) 

^here we use c(-) = cos(-) and s(-) = sin(-) 
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where the eigenvectors are 



V2 



1 

±i sm 7 




(27) 



zti cos 7 

^ It is also possible to represent the whole GRK algorithm as a matrix 

Ggrk ^GiG^'Gi' = (28) 



( ' 



V -1 



3K-4 



K 



It has the form, 



where a and b satisfy - 
element of the group 0(3) 



3K-i 
K(K-l) 







a 
b 




3K-4 

1 
2 



3K-4 
K 



(29) 



(30) 



b^ = 1 which shows that the GRK matrix is an 
9. 



3.1 Reformulation in terms of 5*0(3) Group 

We see in (|26ll that the matrix corresponding to the operator Gi has a reflec- 
tion if ji is odd. This fact makes the analysis of the algorithm difficult, since 
algorithms with even and odd powers of Gi have a different behavior ■ To 
overcome this problem, we can reformulate the algorithm in such a way that 
it will use only even powers of Gi. To do that we have to introduce auxiliary 
search defined by 

Gi = GiG{Gi. (31) 

Now it is necessary to show how the introduction of this new operator is done 
inside the algorithm. We consider a general sequence of Gi and G2 



(u|Gi"Gf-Gf 



'G'^--' Gi^G{'\s{) =0, 



(32) 



We can always make the total number of G\ factors {^jn ~1" jn~2 ~t" ■ ■ ■ J3 ^' .71) 
to be an even number. We can add one extra factor Gi in the beginning using 
the fact that Gi|si) = |si) + 0(1/a/6): 



{u\G{-G'r'G'r'Gi-'^ 



■GfGi'G{'Gi\si) = 0, 



(33) 



We can now consider only sequences with even total number of Gi factors. 
Individual powers of Gi can still be odd. Each odd power of Gi we represent as 
even multiplied by on Gi factor. Since the total number of Gi factors is even. 
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single Gi can only occurs in pairs. This means that in the string of operators 
we can chose single Gi's such that we have 

(7/|---GiG^^---G^^Gi---|si) =0, (34) 

where between the two factors of Gj^ and we have only even powers of 
Gi's. Now one uses the definition of Gq given in Eq. 13 II and rewrites Eq. 1341 as 

H-..G^-.-G^.-.|si> =0. 

Here we used G? = 1 + 0(1/Vfe). 

A general algorithm will now be a sequence of the three operators global 
search Gi raised in even powers, local search G2 and auxiliary search Ga'- 

{u\G{-G'r'&a-' ■ ■ ■ &^&:-G{' = 0. (35) 

Using the matrix description given above we can calculate explicitly a power 
of auxiliary search Ga ■ 

G-J- = 

^ a 

c(2j,02) -c(27)s(2j,02) s(27)s(2j„02) 

c(27)s(2j,02) s2(27) + c2(27)c(2j,02) s(27)c(27)[l - c(2j,02)] 

-s(27)s(2j„02) s(27)c(27)[l-c(2j,02)] 0^(27) + s2(27)c(2j,02) 

(36) 

Its spectrum is 

Gi^ \4) = exp(±2z02j2)l4>, I"2> = l"2> (37) 
Where the eigenvectors can be written as 

\4) = ^i T»cos(27) 1, \u',)=\ sin(27) 1. (38) 



V2 



±i sin(27) 




To summaries in this section: we eliminated reflection by introducing a new 
element Ga. In new formulation of the GRK algorithm (see H35|) only even 
powers of Gi appears. These means that in formulae H26(l and H25|l ji can be 
replaced by even number, so (—1)-'^ = 1. 



4 Lie Algebra Relations 

The introduction of a third operator Ga simplifies the analysis of general al- 
gorithms, since now we do not have to take into account reflections. Now the 
algorithm consists of a sequence of even powers of Gi and integer powers of 
G2 and Ga- Each of these operators [searches ] is an element of SO{3). This 
is a simplification, since now we are dealing with connected component to the 
identity element. But Ga is dependent on the other two operators H31(l . 
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Any element of the SU{2) group can be written in terms of rotations around 
two linearly independent unite vectors (see the Appendix). This is also true for 
5*0(3) group. To find general relations among our three operators it is useful 
to first see what is the Lie Algebra relation. 

Using the matrix form of given by H25|l we can compute its expression 
for small powers: 



sin^ (7) sin (7) cos (7) 

= /+4^?2jTGl = ( - sin^ (7) ) (39) 

- sin 7 cos 7 



The same calculation can be done with G 
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Gi^I + 292 jTc, - / + 292J -1 0, (40) 



and Ga 




— cos (27) sin (27) 
Gi^I + 292jTG^^I + 202j \ cos (27) |, (41) 



note that we used the relation 61 = sin(7) 02 to simplify the above equations. 
The relation follows from the definition of the angles in the limit of 6 ^ 00. 

, T(32 and Tq^ are as elements of Lie Algebra generators corresponding 
to our searches. Using their matrix expressions we see that they satisfy the 
linear relation 

Tg„ + Tg, - 2Tg, = (42) 

which explicitly shows that the Lie algebra elements describing global, local and 
auxiliary searches are linearly dependent. In the next section we shall rise this 
relation into the group, see H45|l . 



5 SU{2) Formulation 

Let us formulate partial search in terms of su{2) algebra and later SU{2) group. 
The transition to SU{2) group makes the manipulation of the group elements 
algebraically easier. 

From Eqs. 1391 1401 and 1411 we see that the Lie algebra generators are linear 
combinations of standard generators and Ty of the so(3) Lie algebra, see (I44f) . 
Any three dimensional vector (|18|l can be mapped to two dimensional matrices 

V V — acTx + bay + cGz , 

using Pauli matrices 

= ( J I), - ( • 7 ) ' '^^ = ( -1 ) ■ (43) 
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We can map so(3) algebra to su{2) algebra by replacing the generators T^, 
Ty and 

/000\ /00-l\ /010\ 

= 1 , Ty = [ , = -1 , 

\0-10/ \100/ \000/ 

(44) 

by |(Ta;, |ctj, and ^a^ correspondingly,so 

R ^ e"-^ ^ u ^ ei" " , 

and rotation act on vectors, 

RV -> uvu-^ 
The actual SU{2) to 0(3) correspondence is given by 

u(r • (t)u^^ = (Rr) ■ a 

Here r is a vector with real components. Elements of partial search can be 
mapped to SU{2) group in the following way: 



cos(ji6li) + i sin(7) sin(ji6'i) - cos(7) sin(ji6li) 

cos(7) sin(ji6'i) cos(ji6'i) - i sin{j) sin(ji6'i) 

r<j2 _ ( e'^'^'^ 

cosUadi) - i cos(27) sin(ja6'2) - sin(27) sin(ja6'2) 

sin(27) sin(ja6'2) cos{ja92) + i cos(27) sin(ja^'2) 

So we mapped partial search in SU{2) group. The elements Ui, U2 and 
of SU{2) group are dependent. Using the Appendix we found algebraic relation 
between group elements describing global, local and auxiliary searches: 

J1„-J2„J1 



uru„-ui' =u^2', (45) 



here sin(7) tan(j2^'2) = tan(ji0i). Corresponding Lie algebraic relation is 142|l . 



6 Conclusion 

In this paper we formulated the partial search algorithms in terms of group 
theory. We think that it will be useful for proof of optimality of GRK algorithm 
in wide class of partial search algorithms. 
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8 Appendix 



Arbitrary element of SU (2) group can be written using only rotation around 
two different axis, see page 176 of the book [14 . We shall represent rotation 
axis by a unit vector. If we define RaW as a rotation around the unit vector ft 
by an angle A: 

-Rn(A) = exp{-i-(n • a)} 

A rotation around any axis by any angle (arbitrary element of SU (2) group R) 
can be represented as sequential rotations around two fixed axis n and m: 

R^Rn{\)R^^ie)Rnh), 

Here n and m are two linearly independent unit vectors and (A, 0, 7) are three 
real numbers (angles). 
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